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Abstract.  The  asymptotic  behavior  of  the  cylindrical  Couette  flow  problem  for  a  rarefied  rotating  gas  with 
evaporation  and  condensation  is  studied  when  the  state  of  equilibrium  is  perturbed  by  the  following  small 
thermodynamic  forces:  (i)  the  pressure  difference  of  the  evaporating  gas;  (ii)  the  angular  velocity  difference 
of  the  cylinders;  and  (iii)  the  temperature  difference  of  the  cylinders.  The  problem  is  solved  by  using  the 
hydrodynamic  equations  that  follow  from  the  balance  equations  of  mass,  momentum  and  energy  of  a  viscous 
heat  conducting  rarefied  gas.  The  hydrodynamic  equations  are  solved  analytically  by  considering  slip  and 
jump  boundary  conditions.  The  fields  of  density,  velocity,  temperature,  heat  flux  vector  and  viscous  stress 
tensor  are  calculated  as  functions  of  the  Knudsen  number  for  each  thermodynamic  force  and  for  different 
values  of  the  angular  velocity.  The  asymptotic  behavior  of  these  fields  are  compared  with  those  obtained 
from  the  kinetic  equation. 


I  INTRODUCTION 

One  of  the  basic  principles  of  continuum  mechanics  is  the  principle  of  material  frame  indifference  which  states 
that  the  constitutive  equations  must  be  the  same  in  inertial  as  well  as  in  non-inertial  systems  of  reference.  This 
principle  has  no  support  within  the  framework  of  kinetic  theory  of  gases  as  was  pointed  out  by  several  authors, 
among  others  we  cite  Muller  [1],  Soderholm  [2],  Sharipov  and  Kremer  [3]  and  Biscari  and  Cercignani  [4]. 

A  very  simple  non-inertial  system  is  represented  by  the  cylindrical  Couette  flow  where  a  gas  is  confined 
between  two  rotating  cylinders.  The  problems  where  the  surfaces  of  the  cylinders  are  rotating  with  different 
angular  velocities  and  are  at  different  temperatures  were  analyzed  in  the  papers  [3].  It  was  found  that  the 
anisotropy  created  by  the  rotation  of  the  gas  induces  the  following  effects:  i)  a  radial  temperature  gradient 
causes  both  radial  and  tangential  heat  fluxes;  ii)  the  diagonal  components  of  the  viscous  stress  tensor  are  not 
zero  but  depend  on  the  rotation  frequency.  These  phenomena  are  observed  not  only  in  weak  but  also  in  strong 
non-equilibrium  Couette  flow. 

The  above  mentioned  papers  consider  the  heat  and  momentum  transfer  only.  To  take  into  account  the  mass 
transfer  we  have  to  assume  that  the  coaxial  cylinders  can  absorb  and  emit  particles  of  the  fluid.  One  of  the 
mechanisms  of  the  mass  exchange  between  the  gas  and  the  cylinders  could  be  the  evaporation  and  condensation 
of  the  particles  on  the  cylinder  walls.  Gas  flows  with  evaporation  and  condensation  on  the  boundary  have  a 
very  complicated  solution  even  in  systems  at  rest,  i.e.  without  rotation.  This  type  flow  was  investigated  by 
many  authors  on  the  basis  of  the  kinetic  theory  and  an  extensive  list  of  the  corresponding  publications  can  be 
found  in  the  paper  [5]. 

The  flow  becomes  more  complex  if  the  cylinders  -  that  confines  a  gas  in  which  processes  of  evaporation  and 
condensation  are  considered  -  rotate.  This  problem  was  solved  recently  by  us  [6]  in  the  whole  range  of  the 
Knudsen  number  by  assuming  that  the  cylinders  have  the  same  temperature  and  rotate  with  the  same  velocity, 
but  they  evaporate  the  gas  with  weakly  different  pressures.  Some  analytical  results  on  the  Couette  flow  with 
evaporation  and  condensation  in  the  continuum  regime  can  be  found  in  the  paper  [7],  while  interesting  results 
on  the  non-linear  Couette  flow  with  evaporation  and  condensation  are  presented  in  the  paper  [8]  where  it  was 
pointed  out  that  a  bifurcation  of  the  flow  is  possible. 

Recently  [9]  the  transport  phenomena  in  rotating  rarefied  gases  that  undergo  evaporation  and  condensation 
on  their  surfaces  were  investigated  by  us  on  the  basis  of  the  kinetic  theory  of  gases.  We  have  considered 
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that  the  state  of  equilibrium  was  perturbed  by  the  following  small  thermodynamic  forces:  (i)  the  pressure 
difference  of  the  evaporate  gas;  (ii)  the  angular  velocity  difference  of  the  cylinders;  and  (iii)  the  temperature 
difference  of  the  cylinders.  For  every  thermodynamic  force  all  transport  phenomena,  of  mass,  momentum  and 
heat  transfer  were  calculated,  i.e.  the  direct  effects  such  as:  (i)  the  mass  transfer  caused  by  the  pressure 
difference,  (ii)  the  momentum  transfer  caused  by  the  velocity  difference,  and  (iii)  the  heat  flux  caused  by  the 
temperature  difference.  Further  the  Onsager-Casimir  reciprocity  relations  for  the  cross  effects  were  verified. 
The  used  kinetic  equation  was  the  Shakov  model  of  the  Boltzmann  equation  and  it  has  been  solved  by  the 
discrete  velocity  method  where  the  discontinuity  of  the  distribution  function  has  been  taken  into  account. 

The  aim  of  this  work  is  to  And  the  asymptotic  behavior  of  the  above  mentioned  phenomena  by  using  the 
hydrodynamic  equations  that  follow  from  the  balance  equations  of  mass,  momentum  and  energy  of  a  viscous 
heat  conducting  rarefied  gas.  The  hydrodynamic  equations  are  solved  analytically  by  considering  slip  and  jump 
boundary  conditions  [10].  The  fields  of  density,  velocity,  temperature,  heat  flux  vector  and  viscous  stress  tensor 
are  calculated  as  functions  of  the  Knudsen  number  for  each  thermodynamic  force  and  for  different  values  of  the 
angular  velocity.  The  asymptotic  behavior  of  these  fields  are  compared  with  those  obtained  from  the  kinetic 
equation  [9]. 


II  HYDRODYNAMIC  SOLUTION 


We  consider  a  rarefled  gas  between  two  coaxial  rotating  cylinders  with  radii  Rq  and  i?i  {Ro  >  Ri).  The 
cylinders  are  rotating  around  the  x-axis  and  they  are  supposed  to  be  so  long  that  end  effects  can  be  neglected. 
Further  we  assume  that  the  equilibrium  state  is  weakly  disturbed  by  the  following  three  factors:  (i)  The  outer 
cylinder  evaporates  the  particles  with  an  equilibrium  pressure  Po{Ro),  while  the  pressure  of  evaporated  particles 
from  the  inner  cylinder  Pi  slightly  differs  from  the  equilibrium  pressure  at  its  surface,  i.e. 

P,  =  PoiRi)  +  AP,  (1) 


The  equilibrium  pressure  is  given  by  the  equation  of  state  Po(r')  =  no{r')kTo  where  r'  =  is  the  radial 

coordinate  and  no(r')  is  the  particle  number  density  between  the  cylinders 


no(r') 


noo(l  -  RyR^)mM^ 

exp[(,5flo-Ro)^]  -  exp[(,dnoJ?i)^] 


exp[(^fior')^], 


(2) 


In  the  above  equation  noo  denotes  the  number  density  when  the  cylinders  are  at  rest  (flo  =  0),  m  is  the  mass 
of  a  gas  particle  and  k  is  the  Boltzmann  constant;  (ii)  The  outer  cylinder  rotates  with  an  angular  velocity  Uq, 
while  the  angular  velocity  of  the  inner  cylinder  fli  slightly  differs  from  Qq,  i.e. 


ni=flQ  +  Afi,  j3Pi|An|<l. 


(3) 


The  smallness  of  AO  means  that  the  difference  between  the  velocities  of  the  surfaces  is  smaller  than  the  sound 
speed;  (iii)  The  outer  cylinder  is  at  the  equilibrium  temperature  Tq,  while  the  temperature  of  the  inner  cylinder 
Ti  slightly  differs  from  Tq,  i.e. 


Ti  =  To  +  AT,  ^  <C  1.  (4) 

Jo 

There  are  two  main  parameters  that  determine  the  solution  of  the  problem:  the  inverse  Knudsen  number 
which  is  the  rarefaction  parameter  S  =  v^Po/(2Aoo)  where  Aqo  is  the  molecular  mean  free  path  at  the  density 
noo  and  at  temperature  To;  and  the  dimensionless  angular  velocity  w  =  /30oPo- 

In  order  to  solve  the  above  mentioned  problem  in  the  hydrodynamic  regime  flow  we  use  the  balance  equations 
of  mass  density  g,  momentum  density  gu'^  and  internal  energy  density  ge  that  read 


dt  dx', 


=  0, 


dgu'i  d  ,  ,  ,  „  dge  d  ,  ,  ,,  ,  du',  „ 

+  ~^ieUiUj+P8ij  -Qij)  +  5i)  -  0-  (5) 


These  balance  equations  are  supplemented  by  the  constitutive  equations  for  the  viscous  stress  crb,  heat  flux 
vector  q[  and  specific  internal  energy  e: 
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,  du'i  du'j  2  du'^ . 

J  %  1 


Qi  =  -«l 


The  above  equations  represent  the  constitutive  equations  of  a  Newtonian  heat  conducting  fluid  without  bulk 
viscosity.  The  coeflicient  of  thermal  conductivity  k  is  given  in  terms  of  the  coefficient  of  shear  viscosity  fx  by 
K  =  15fcju/(4m).  For  a  monatomic  ideal  gas  k  and  fx  are  only  functions  of  the  temperature  T. 

By  taking  into  account  that  all  fields  depend  only  on  the  radial  coordinate  r'  and  that  we  are  looking  for  a 
stationary  solution  of  the  gas  flow,  Eqs  (5)  reduce  to 


r'dr'XT  n 


/  ,dK  < 


dP  4  d  \ ^x  d{r' u'^)!  2  ,  dfx 

n  I  d”  n  I  f  I  I  I  ^T*  t  ; 


dr'  3  dr'  r'  dr' 


/  ^  ^  r  ,3 A  ( < Y 

^  V  ’■  dr'  r'  )  r'2  dr'  T  dr'  \r'  )\' 


-  4^  +4<  +  \r'±  (<']]]  (9) 

2T  ’’dr'  r'dr'V  dr'  ’'dr'  ^  |  3  r'  dr'  r'2  dr'  r'^  dr'  \  r' J  |'  ^  ^ 


3  r'  dr' 


dr'  V  r 


In  the  above  equations  the  mass  density  g  has  been  represented  by  the  equation  of  state  g  =  mP/kT.  Eqs  (7) 
-  (9)  refer  to  a  system  of  four  coupled  differential  equations  for  the  fields  u'r(r'),u'^(r'),  P(r'),  Tir')  which  we 
proceed  to  And  its  solution  for  given  boundary  conditions. 

We  note  first  that  for  a  linearized  problem,  in  which  the  conditions  (1),  (3)  and  (4)  hold,  we  can  write 

P(r')=Po(r')  l  +  ,  <(r')=r'  ,  (10) 

a  J  La. 

u;(r')  =  nor'+r'  ’  T{r')=To  1  +  ,  (11) 

-  a  J  La. 

where  a  =  P,fl,  T,  v^°‘^ ,  ui°‘^ ,  and  are  dimensionless  quantities  and  are  the  thermodynamic  forces 

Xp  =  AP/Po(i?i),  Xn  =  and  Xp  =  AT/Tq. 

Next  we  insert  (10)  and  (11)  into  Eqs  (7)  -  (9)  and  get  the  linearized  field  equations  for  the  dimensionless 
quantities  and  that  read 


dln[norUr“^] 


r  -  4a;u^  +  3no  "  (5  ’ 


dr  r  dr  nnn 


dr  [  dr  15no  ^  ’ 


where  r  =  r'/Po  is  a  dimensionless  quantity.  The  general  solution  of  the  system  (12)  -  (13)  for  the  fields 
and  are 


=  — — noo, 
nor  ^ 


=  2SuD^°‘^  Inr  —  P{“V  +  — — ^  Inr  +  (7^“^  —  , 


=  85w2p(“)  ^g{r)  -  +  5^“^  inr  +  (r^  -  2(7^“V^  -  2C[°‘'>g{r)^  + 


4  nlo  P(“)a;2 
3ng(r)  S 
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where  b[‘^\  and  are  constants  of  integration  and  g{r)  =  Inr  —  r^/2. 

Once  the  fields  ui°‘\r),  u^\r),  r^“^(r),  v^°‘\r)  are  known  we  can  also  obtain  the  non- vanishing  components 
of  the  heat  flux  vector  and  for  the  viscous  stress  tensor  from  (6),  (14)  and  (15): 


3'“^  = 


15  1 

8  6  dr 


2 

3  dr  3  r 


1  ui°‘^ 


1  £)(“) 


noo 


3  dr^  no{r) 


(3  -I-  4w^r)  , 
(16) 


j(“)  = 

w 


2  1 


3  r 


3  dr 


1 


^00 


3  Sr"^  no(r) 


(1  -  2wV)  ,  ^ 


du!^^ 

dr  r 


=  1-2 


B. 


(“) ' 


For  w=0  the  solution  of  the  system  of  differential  equations  (12)  and  (13)  reduces  to 


u[“V  =  D^°‘\ 


=  e[°‘\  -h  Inr  -f  =  G^°‘\ 


(a) 


(17) 


(18) 


where  D^°‘^ ,  e[°‘^  ,  E^^ ,  and  are  constants  of  integration. 

Ill  SLIP  AND  JUMP  BOUNDARY  CONDITIONS 


In  order  to  determine  the  constants  of  integration  of  the  previous  section  one  has  to  know  the  boundary 
conditions  at  the  outer  cylinder  (r  =  1)  and  at  the  inner  cylinder  (r  =  Ri/Ro).  The  slip  and  jump  boundary 
conditions  adopted  here  that  correspond  to  a  perturbation  of  the  pressure,  temperature  and  angular  velocity 
on  the  outer  cylinder  are: 


t)(“)  = 


S  no(l)  dr 


r(“)=a5uW-^^^-2 


(«) 


CZf  ?^00 

S  no(l)  dr  S  no(l)  dr 


(a)  “9  »^00  f  du\f^  A 

while  the  corresponding  slip  and  jump  boundary  conditions  on  the  inner  cylinder  read 

,(«) 


04  ?^00  (a) 

6  no(l)  "  ’ 

(19) 

as  »^00  (a) 

s  noil)  ^  ’ 

(20) 

(21) 

;(“)  =  d 


Pa  OCiilj. 


(a)  ^2  Tloo 

r  "1  r.  /  —  t  . 


dT^°‘'>  ^2“^  ”00  du^r 


=  dra  -  Q:6U^“^  -I-  ^ 


6  no(i7i/i?o)  dr  6  no(J?i/i?o)  dr 

(«) 


0:7  npo  dur 

4“  2  - 


04  noo  (a) 
S  no(i?i/iJo)  ’ 


”oo_  ^(«) 


S  no{Ri/Ro)  dr  6  no(Ri/Ro)  dr  S  no{Ri/Ro)  ’’  ’ 


(22) 


(23) 


09  *^00 


du, 


=  4o  -I- 

^  (5  no{Ri/Ro)  \  dr 


(a) 

V 


,(«)' 

l(p 

r 


(24) 


In  the  above  equations  Spa,  Spa  and  fca  are  Kronecker  symbols  and  a  =  P,T,  0. 
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The  terms  in  the  jump  of  the  pressure  as  stated  in  the  boundary  conditions  above  have  the  following  meaning: 
(a)  is  due  to  the  evaporation  and  condensation  at  the  solid  boundary;  (b)  {a2/d){dT^°‘^  /dr)  is  due  to 

a  normal  gradient  of  temperature  at  the  solid  boundary;  (c)  (2a3/S)(dui°‘^ /dr)  is  due  to  a  normal  gradient  of 
the  normal  component  of  the  velocity  at  the  solid  boimdaxy;  (d)  {a4/S)ui°‘^  and  (2a4/5)ur“^  are  due  to  the 
curvature  of  the  solid  surfaces.  The  modulus  of  the  mean  curvature,  as  defined  in  the  work  by  Sone  [10],  is  equal 
to  one  for  the  internal  cylinder  and  1/2  for  the  external  cylinder.  The  terms  in  the  jump  of  the  temperature 
have  the  same  meaning  as  those  that  appear  in  the  jump  of  the  pressure.  The  jump  in  the  tangential  velocity 
is  due  to  the  normal  gradient  of  the  tangential  velocity.  The  coefficients  ai  through  0:9  are  calculated  from  a 
kinetic  equation  and  the  values  for  these  coefficients,  based  on  the  BGK  model,  are  presented  by  Sone  [10]. 
Since  we  have  applied  the  S  model  to  solve  the  present  problem  in  the  kinetic  regime  [9],  we  have  to  use  the 
coefficients  based  on  the  S  model  too.  It  was  verified  that  the  coefficients  and  og  based  on  the  S  model 
are  3/2  times  greater  than  those  based  on  the  BGK  model,  while  the  other  coefficients  are  exactly  the  same. 
Thus,  the  numerical  data  presented  by  Sone  [10]  can  be  used  here  if  the  two  above  mentioned  coefficients 
are  multiplied  by  3/2.  We  have:  =  2.13204,  0:2  =  0.83766,  0:3  =  0.82085,  04  =  0.38057,  og  =  0.44675, 

ag  =  1.95408,  07  =  0.33034,  og  =  0.13157  and  =  1.01619. 

Once  the  boundary  conditions  are  known  it  is  easy  to  obtain  the  constants  of  integration  which  appear  in 
the  hydrodynamic  solutions  (14),  (15)  and  (18).  Hence  the  fields  of  density,  velocity,  temperature,  heat  flux 
vector  and  stress  tensor  can  be  determined  as  functions  of  the  rarefied  parameter  and  of  the  angular  velocity. 


IV  RESULTS 

For  the  fiow  caused  by  a  pressure  difference  we  have  plotted  in  Figures  1  to  3  the  kinetic  and  the  hydrodynamic 
solutions  of  the  fields  and  as  functions  of  the  rarefied  parameter  S  for  different  values  of  the 

angular  velocity  w.  We  note  that  a  better  agreement  of  the  hydrodynamic  and  kinetic  solutions  for  and 
qr^^  is  attained  for  low  angular  velocities.  The  difference  in  the  hydrodynamic  and  kinetic  solutions  for 
is  due  to  the  fact  that  it  is  small  as  any  cross-effect.  The  kinetic  and  the  hydrodynamic  solutions  of  the  fields 
and  (T-pip  for  a  fiow  caused  by  the  velocity  difference  are  plotted  in  Figures  4  to  6  as  functions  of  the 
rarefied  parameter  5  and  of  the  angular  velocity  oj.  While  there  is  a  good  agreement  of  the  hydrodynamic  and 
kinetic  solutions  for  ,  the  difference  between  the  solutions  for  and  qi^^  is  due  to  the  fact  that  both 
represent  fields  of  cross-effects,  which  are  usually  small.  The  fields  u^\  qi^^  and  which  are  solutions  of 
the  hydrodynamic  and  kinetic  equations  are  plotted  in  Figures  7  to  9  for  a  flow  caused  by  the  temperature 
difference  as  functions  of  the  rarefied  parameter  6  and  of  the  angular  velocity  uj.  There  is  a  very  good  agreement 
of  the  hydrodynamic  and  kinetic  solutions  for  ai-ip'  and  q)-  ’  and  we  note  that  for  high  rarefaction  parameter 
the  heat  flux  q)  ’  does  not  depend  on  w  and  all  curves  coincide.  The  field  u'r  '  corresponds  also  to  a  cross-effect 
but  here  the  agreement  between  the  solutions  is  better  that  the  two  cases  analyzed  above. 
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FIGURE  1.  Ur  ^  as  function  of  S  for  different  values  of  oj 
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FIGURE  2.  as  function  of  6  for  different  values  of  w, 
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FIGURE  3.  as  function  of  5  for  different  values  of  w 
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